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We provide the theoretical basis for understanding the phenomenon in 
which an ultra cold atom incident on a possibly warm target will not stick, 
even in the large n limit where n is the number of internal degrees of freedom 
of the target. Our treatment is non-perturbative in which the full many-body 
problem is viewed as a scattering event purely within the context of scatter- 
ing theory. The question of sticking is then simply and naturally identified 
with the formation of a long lived resonance. One crucial physical insight 
that emerges is that the many internal degrees of freedom serve to decohere 
the incident one body wavefunction, thus upsetting the delicate interference 
process necessary to form a resonance in the first place. This is the physical 
reason for not sticking. 



I. INTRODUCTION 

The problem of low energy sticking to surfaces has attracted much attention over the 
years [0-^. The controversial question has been the ultralow energy limit of the incoming 
species, for either warm or cold surfaces. A battle has ensued between two countervailing 
effects, which we will call classical sticking and quantum reflection. The concept of quantum 
reflection is intimately tied into threshold laws, and was recognized in the 1930's by Lennard- 
Jones [|l[]. Essentially, flux is reflected from a purely attractive potential with a probability 
which goes as 1 — a^/e, as e — > 0, where a is a constant and e is the translational energy 
of the particle incident on the surface. Classically the transmission probability is unity. 
Reflection at long range prevents inelastic processes from occurring, but if the incoming 
particle should penetrate into the strongly attractive region, the ensuing acceleration and 
hard collision with the repulsive short range part of the potential leads to a high probability 
of inelastic processes and sticking. 

The blame for the quantum reflection can be laid at the feet of the WKB approximation, 
which breaks down in the long range attractive part of the potential at low energy. Very far 
out, the WKB is good even for low energy, because the potential is so nearly flat. Close in, 
the kinetic energy is high, because of the attractive potential, even if the asymptotic energy 
is very low, and again WKB is accurate. But in between there is a breakdown, which has 
been recognized and exploited by several groups P-pl- We show in the paper foiling this 



one that the breakdown occurs in a region around \V\ ~ e; i.e. aproximately where the 
kinetic and potential energies are equal. 

It would seem that quantum reflection would settle the issues of sticking, since if the 
particle doesn't make it in close to the surface there is no sticking. (Fig [I]) There is one 
caveat, however, which must be considered: quantum reflection can be defeated by the 
existence of a resonance in the internal region, i.e. a threshold resonance. (Fig |2]) 

The situation is very analogous to a high Q Fabry-Perot cavity, where using nearly 100% 
reflective, parallel mirrors gives near 100% reflection except at very specific wavelengths. 
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At these specific energies a resonace buildup occurs in the interior of the cavity, permiting 
near 100% transmission. Such resonances are rare in a one dimensional world, but the huge 
number of degrees of freedom in a macroscopic solid particle makes resonance ubiquitous. 
Indeed, the act of colliding with the surface, creating a phonon and dropping into a local 
bound state of the attractive potential describes a Feshbach resonance. Thus, the resonances 
are just the sticking we are investigating, and we must not treat them lightly! Perhaps it is 
not obvious after all whether sticking occurs. 

After the considerable burst of activity surrounding the sticking issue on the surface 
of liquid Helium [T2|jr3|] , and after a very well executed theoretical study by Clougherty 



and Kohn [Q, the controversy has settled down, and the common wisdom has grown that 
sticking does not occur at sufficiently low energy. While we agree with this conclusion, 
we believe the theoretical foundation for it is not complete, nor stated in a wide enough 
domain of physical situations. For example, Ref. J|| treats only a harmonic slab with one 
or two phonon excitation. It is not clear whether the results apply to a warm surface. On 
the experimental side, even though quantum reflection was observed from a liquid Helium 
surface, that surface has a very low density of available states (essentially only the ripplons) 
which could be a special case with respect to sticking. Thus, the need for more rigorous 
and clear proof of non-sticking in general circumstances is evident. This paper gives such 
an analysis. In a following paper, application is made to specific atom-surface and slab 
combinations, and the rollover to the sticking regime as energy is increased (which can be 
treated essentially analytically) is given. 

The strategy we use puts a very general and exact scattering formalism to work, providing 
a template into which to insert the properties of our target and scatterer. Then very general 
results emerge, such as the non-sticking theorem at zero energy. The usual procedure of 
defining model potentials and considering one phonon processes etc. is not necessary. All 
such model potentials and Hamiltonians wind up as parameters in the R-matrix formalism. 
The details of a particular potential are of course important for quantitative results, but the 
range of possible results can be much more easily examined by inserting various parameters 
into the R-matrix formalism. All the possible choices of R-matrix parameters give the correct 
threshold laws. Certain trends are built into the R-matrix formalism which are essentially 
independent of the details of the potentials. 

Before commencing with the R matrix treatment, we briefly consider the problem per- 
turbatively in order to better elucidate the role played by quantum reflection. We emphasize 
that none of the perturbation section is actually necessary for our final conclusions. 

In a perturbative treatment for our slab geometry, quantum reflection simply results in 
the entrance channels' wave function (at threshold) having its amplitude in the interaction 
region go to zero as k e ~ y/e when normalized to have a fixed incoming flux. (k e is the 
magnitude \k e \ of the incident wavevector of the incoming atom). The inelastic transition 
probabilities are proportional to the potential weighted overlap of the channel wavefunctions 
and this immediately leads to the conclusion that the inelastic probability itself vanishes as 
k e ~ \/e- As mentioned, this conclusion is shown to rigorously remain true using the R 
matrix. We show in this paper that in spite of the inherently many-body nature of the 
problem, in the ultra-cold limit we can correctly obtain the long-range form of the entrance 
channel's wavefunction by solving for the one- dimensional motion in the long-range surface- 
atom attraction (i.e. the diagonal element of the many-channel potential matrix). This 
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FIG. 1. The stationary state one body wavefunction of the incident atom moving in the 
y-independent mean potential felt by it. The amplitude inside the interaction region is supressed 
by k e ~ sfi. This is tantamount to the reflection of the atom. 




target absorbs 
energy 



FIG. 2. A schematic view of a Feshbach resonance wherein the incident atom forms a long lived 
quasi-bound state with the target. The many body wavefunction in this situation (not shown) has 
a large amplitude in the 'interior' region near the slab. 

allows quantitative predictions of the sticking probability, which we do in the following 
paper. There, we further exploit the perturbative point of view together with an analysis of 
WKB to predict a 'post-threshold' behavior as quantum reflection abates, when the incoming 
energy is increased. 

II. GEOMETRY AND NOTATION 

The incident atom is treated as a point particle at position (x,y). To keep the notation 
simple we leave out the z-coordinate and confine our discussion to two spatial dimensions. 
Thus a cross-section will have dimensions of length etc. It will be quite obvious how and 
where z may be inserted in all that follows. Let u represent all the bound degrees of 
freedom of the scattering target, which we take to be a slab of crystalline or amorphous 
material. Let fl c (u), c — 1, 2, • • •, be the manybody target wave functions in the absence of 
interactions with the incident particle, and having energy £' arget . These are normalized as 
/ail u \VL c {u)\ 2 = 1. x is the distance of the scatterer (atom) from the face of the slab which 
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is approximately (because the wall is rough) along the line x — 0. The internal constituents 
of the slab lie to the left of x = and the scatterer is incident from the right with kinetic 
energy e = h 2 k 2 /2m. The total energy E of the system is 

E = e + Ef rget (1) 

where c = e is the index of the 'entrance channel' i.e. the initial internal state of the slab 
before the collision is Q e (u). Notice that we say nothing about the value of E^ vget itself. 

— * 

In particular the slab need not be cold. k c is the magnitude of the wave vector k c of the 
particle when it leaves the target in the state Q c ( u ) after the collision. Our interest focusses 
on k e — > 0. k e is the magnitude of the wavevector of the incoming particle. For the open 
channels c = 1 , • • • n (this defines n) for which E > El arget 

^ f"<*-5' ga) (c<»); (2) 

whereas for the closed channels (c > n), E < E^ 561 and 

, • /2m(^ arget - E) . . 

k c = «y ^2 = lK c {On). (3) 

k c > 0. We will use {k cx ,k cy ) as the x, y components of k c . Let Ui nt {x,y,u) = 
{2m / h 2 )Vi nt {x , y, u), where Vi nt {x,y,u) describes quite generally the interaction potential 
between the incident atom and all the internal degrees of freedom of the slab. For simplicity 
we assume for the moment that there is no interaction between slab and atom for x > a. 



III. PRELIMINARIES: PERTURBATION 

As stated above, we excercise the perturbative treatment for insight only; our final 
conclusions are based on nonperturbative arguments. 

We treat the interaction U int {x,y,u) between slab and atom by separating out a 'mean' 
potential felt by the atom that is independent of y and u; call it U^{x). The remainder 
U^{x,y,u) = U- mt {x,y,u) — U^{x) is treated as a perturbation. 

Now the incident beam is scattered by the entire length (say from y = —L to L = 2L) of 
wall which it illuminates. If all measurements are made close to the wall so that its length 
2L is the largest scale in the problem, then it is appropriate to speak of a cross-section per 
unit length of walk a dimensionless probability. More specifically, we will assume that the 
matrix elements U C J (x, y) = / du n*{u)U^{x, y, u)Q c >{u) of the perturbation U^{x, y, u) 



all u 



in the Q c (u) basis are given by the simple form U c \){x ) y) = U c \J{x)f{y) for y e [— L, L] 
and elsewhere. f{y) is a random persistent (does not die to as \L\ — > oo) function 
that models the random roughness of the slab and is characterized by its so-called spectral 
density function S, a smooth positive-valued non-random function, such that 



Li 

J dye ik yf{y) 

-L 



= 2LS{k) \/k 



(4) 
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as L — > oo. 

Now, applying either time-independent perturbation (equivalently the Born approxima- 
tion for this geometry) or time-dependent perturbation theory via the Golden Rule, gives 
that the cross-section per unit length of wall for inelastic scattering to a final channel c is 

P%-e(0) = 2 -r(l dx'<t>{x'\ k cx )U£\x')<P(x'; k ex )) S(k cy - k ey ) (5) 



-oo 



where <f>(x\ k x ) is the solution of the o.d.e. 

(J^-U^(x) + k 2 }j ( p(x;k x ) = (6) 

which is regular or goes to zero as x — > — oo inside the slab and is normalized as 

4>(x; k x ) ~ sm(k x x + S) asx — > oo (7) 

Accepting for the moment that as k e — > the amplitude of (p(x; k ex ) in the internal region 
x < a goes to zero as k e ~ y/e, then the square of the overlap integral in Eq. behaves 
as fcg, because by our proposition the amplitude of 4>{x'\ k ex ) ~ k ex ~ k e . Together with the 
l/k e prefactor we get an overall behavior of k e for the inelastic probability as claimed. 
To show that indeed as k e — > the amplitude of 4>(x; k ex ) in the internal region x < a goes to 
zero as k e ~ \/e, we temporarily disregard the required normalization of <J)(x; k x ) of Eq. ([7[) 
and fix its initial conditions (slope and value) at some point inside the interaction region 
x < a such that the regularity condition is ensured. We then integrate out to x — a. Let 
us denote this unnormalized solution with a prime, as (p'(x; k x ). The point is for k x varying 
near 0, both v (the value) and s(the slope) that the solution emerges with at x — a, are 
independent of k x and in fact the interior solution thus obtained is itself independent of k x . 
This is because the local wave vector k(x) = \j2m(e — U(x))/h 2 essentially stays the same 
function of x for all e near 0. Therefore for x > a 4>(x; k x ) continues onto 



v 



cos[k x (x — a)] + — sin[k x (x — a)} x > a (8) 
k x 



This is a phase-shifted sine wave of amplitude ~ l/k x . We must enforce the normalization 
of Eq. (|7|) and get <p(x; k x ) ~ k x cj)'(x; k x ). As a result, the interior solution gets multiplied by 
k x and we thereby have our result. <p(x; k x ) is the solution of a one-dimensional Schrodinger 
equation for the incoming particle in the one- dimensional long-range potential created by 
the slab. The suppression of its amplitude by \fl near the slab is due to the reflection 
it suffers where the interaction turns on. Within the perturbative set-up the non-sticking 
conclusion is then already foregone 0. 

The problem is whether we can really accept this verdict of the one-dimensional un- 
perturbed solution, when in fact we know that the turning on of the perturbation (many 
body interactions) causes a multitude of resonances to be created, internal resonances being 
exactly the situation in which the Proposition above is known to badly fail. It appears that 
the perturbation is in no sense a small physical effect. Therefore a nonpeturbative approach 
is needed. Here we use R-matrix theory in its general form to accomplish the task. 
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IV. S-MATRIX AND R-MATRIX 



One point that the preceding section has made clear is that it is the energies (both initial 
and final) in the x-direction, perpendicular to the slab that are most relevant. In fact as 
regards the final form of our answers the motion of the y degree of freedom may as well 
have been the motion of another internal degree of freedom of the slab. In other words, 
mathematically speaking, the y degree of freedom may be subsumed by incorporating it 
as just another u. For example, we may imagine the incident atom being confined in the 
y-direction by the walls of a wave-guide at y — — LandL that is large enough so that it could 
not possibly change the physics of sticking. Then we quite rigorously have a bound internal 
state of the form 

Q Cjn (y,u) = Q c (u) sin -^- (9) 

x is now the only scattering degree of freedom. There will be no necessity in carrying along 
the extra index n and variable y as in Eq. (H), and we will simply continue to write fl c (u) 
instead. Thus with this understanding, the problem is essentially one-dimensional in the 
scattering degree of freedom. 

We proceed to derive the expression for the S matrix in terms of the so-called R matrix, 
and derive the structure of the R matrix. For simplicity we continue to assume for the 
moment that there is no interaction for x > a. Then for x > a, the scattering wavefunction 
of the interacting system corresponding to the scattering particle coming in on one entrance 
channel, say c = e, with energy e = h 2 k 2 /(2m) is 

oo / e ~ik e x e ik c x \ 
V>0, U) = \ n- $ce pr S ce ) ^c{u) X> d (10) 



c=l \ V ,v e V ,v c / 

where the sum must include all channels, even though the open channels are finite in number. 
The factors of k~ x l 2 in Eq. fllQf) mean that the flux in each channel is proportional only to 
the square of the coefficient and hence ensure the unitarity of S. With this convention, the 
open-open part of the S- matrix — the nxn sub matrix S cc > with c, d — 1, 2, . . . , n — is unitary. 
\fk~ c = e m l A ^/~K~ c may be arbitrarily chosen since it cannot affect the open-open part of S. 

S is found in analogy to the one- dimensional case by introducing the matrix version of 
the inverse logarithmic derivative at x = a called R(-E) the Wigner R-matrix defined by 

v = R(E)s (11) 

where the components of v and s are the expansion coefficients of ip(x = a,u) and <tt^!t±m 
respectively in the Vt c {u) basis. Supposing 8 ^^~ a ' M - > to be known, we will (like in electro- 
statics) use the Neumann Green's function Gn{x,u]x',u') to construct ip(x,u) everywhere 
in the interior x < a. ip(x, u) satisfies the full Schrodinger equation with energy E. We need 
Xx(x, u) A = 1, 2, • • • , the normalized eigenfunctions of the full Schrodinger equation in the 
interior x < a with energies E\, satisfying Neumann boundary conditions 8x< -g~ ct ' n ' ) = 0. So 

V 2 + V int (x,u)-E)^(x,u) = (12) 



2m 
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h 2 \ 

V 2 + V int ( Xj u) - E x ) X x(x,u) = (13) 



2m 
2m 



V 2 + V int (x, u)-E) G N (x, u; x', u') = 8{x - x')8{u - v!) (14) 



where V 2 = Jj + and 

gg^^V) =0 anrf ax(x = a, M ) =0 (15) 

ax ax 
=>G N (x,u;x,u) = 2^ (16) 

A=l 

Gn is symmetric in the primed and unprimed variables. By Stokes' Theorem, 

(~h 2 /2m) J dx' J du' (0i V ,2 2 - 2 V%) = (~h 2 /2m) J du' (0i V^,0 2 - <faVM 

x'<a all u' x'=a, all u' 

where V^(-) = £'(•) ■ V with <pi = if)(x',u') and 2 = Gn(x,u;x',u') gives 

,/ x ^ f j , ^ i i r >dif)(x' = a,u') . . 

w[x,u) = / du Gn(x } u;x } u) x<a (18) 

2m J ox' 

all u' 

Put x = a and it is deduced using Eqs. (|TTD and together that 



a=i E * ~ E 



(17) 



RcAE) = Y,^^ (19) 
where 7a c = J ofo xaK u)Q c (u). 



all u 



A. The S matrix 

Now shifting attention to the outside (x > a), we see that we can compute both V^V( a > u ) 
and ^(a, u) on the surface a; = a using the asymptotic form of Eq. (|T0| ) which automatically 
gives these expanded in the £l c (u) basis. Writing the matrix Eq. ( O ) is now simple. It 
is best to do it all in matrix notation, and thus be able to treat all possible independent 
asymptotic boundary conditions simultaneously. 

Let e lkx , \/~k and 1 / \f~k be diagonal matrices with diagonal elements e lkcX } \/k^ and 1 / \fk c . 
Then Eq. (|TT| ) reads 

p —ika p ika / p —ika p ika \ 

S = iRh ( ——j= —t=S ] . (20) 



Each column c = 1, ... , n of the matrix equation above is just Eq. (ITT]) for the solution 
corresponding to an incoming wave only in channel c (For c> n the wavefunctions blow up 
as x — > oo). Remembering that non-diagonal matrices don't commute, we solve for S to get 
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S = e~ ika Vk — -(1 + iRk)^=e~ ika (21) 



or, with some simple matrix manipulation, 



S = e~ lka — -={l + iVkRVk)e~ ika . (22) 

1 2 v /lIvv k 



V. S MATRIX NEAR A RESONANCE 

As discussed in the introduction, the resonances are a key to the sticking issue. Sticking 
is essentially a long lived Feshbach resonance in which energy has been supplied to surface 
and bulk degrees of freedom, temporarily dropping the scattering particle into a bound 
state of the attractive potential. Thus we must study resonances in various circumstances 
in the low incident translational energy regime. We derive the approximation for S(E) near 
E = Eq, a resonant energy of the compound system. Eq is the total energy of the joined 
(resonant) system. Within the R- matrix approach, the Xx{ x -> u ) °f section |V| are bound, 
compound states with Neuman boundary conditions at x — a. .R-matrix theory properly 
couples these bound state to the continuum, but some of the eigenstates are nonetheless 



weakly coupled to the continuum, as evidenced by small values of the 7a c 's of section [TV]; 
these are the measure of the strength of the continuum couplings. While every one of the 
R— matrix bound states will result in a pole E\ in the R matrix expansion, only the weakly 
coupled ones are the true long lived Feshbach resonances of physical interest. It is also 
helpful to know that the values of these 'truly' resonant poles at E\ are the most stable to 
changes in the position x = a of the box. This in fact provides one unambiguous way to 
identify them. Our purpose here is to derive the resonant approximation to the S matrix in 
the vicinity of one of these Feshbach resonances. We do so using the form of the R-matrix 
in Eq. dT9|) . Note that the energy density p(E) = 1/D(E) of these Feshbach resonances will 
be large because of the large number of degrees of freedom of the target. D(E) is the level 
spacing of the quasibound, resonant states. 



A. Isolated Resonance 

As mentioned, the point of view we will take is to identify a resonant energy with a 
particular pole E\ in the R matrix expansion of Eq. ([19|) . Those E\ corresponding to 
resonances are a subsequence of the E\ appearing in the expansion in Eq. flT5|). For E near 
a well isolated resonance at E\ we separate the sum-over-poles expansion of the R-matrix 
into a single matrix term having elements ^|^rff, plus a sum over all the remaining terms, 
call it N. If the energy interval between E\ and all the other poles is large compared to the 
open-open residue at E\ then we may expect that the nx n open-open block of iV will have 
all its elements to be small. Then rewriting the inverse in Eq. fl2"2"|) 
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where M = \fkN\fk and V cc i = (v / ^c7ac)(v^V7a c / )) an d setting M = allows us to simplify 
the central term in Eq. (1221) exactly. (We will return to the case M ^0.) 



1 



1 

1 + 
1 + 

1 + 



1 



Ex-E 
1 



1 + iv^fcRVfc) 
2iVkRVk 
2/ — — - (withM = 0) 



i\/kTlVk 
1 _. V 



E-iV 
1 



E x -E 
-2iV 



E x -E- i(T x /2 + tAE) 



2iVk 



where we used 



V 2 



lln^n) + (7A(n+l) K "+l H )) 

■) + *(7A(n+l) K n+l + •■•)) 



2 



+ ••• + 



iAE\ V 



V 
V 



to get the identities 



[Ex-E- iV]V = [E X -E- i(T x /2 + tAE)]V 
1 „ 1 



V 



Ex-E- i(T x /2 + iAE) E x - E - iV 



:V 



(24) 
(25) 
(26) 

(27) 
(28) 

(29) 
(30) 

(31) 

(32) 
(33) 



Also define (Txc/2) 1 ^ 2 = jxcVkc, c — 1, 2, • • • , n. This defines the sign of the square-root on 
the lhs. to be the sign of 7a c and allows the convenience of expressing things in terms of the 
I\ c 's and their square-roots, and not having to use the 7a c 's themselves. Thus we arrive at 



S r .r' 



-ik c a 



T l/2 r l/2 
LL Ac 1 Ac' 



-ik n ia 



Et ] -E- iTx/2 



(34) 



(r) 

where E x = E x + AEx, for the n x n open-open unitary block of S in the neighbourhood 
of a single isolated resonance after neglecting the contribution of the background matrix M. 
For us the essential point is that 



r Ac = 2 k c (E) 7 2 Xc 



(35) 



that the partial widths Tx c depend on the energy E, through the kinematic factor k c (E). 
Mostly this energy dependence is small and irrelevant except where the /c c 's and hence T^'s 
are varying near 0. These are the partial widths of the open channels near threshold. Hence 
I'S'cel 2 (c 7^ e) an inelastic probability behaves like k e ~ y/e when the entrance channel is at 
threshold. Including the background term (M ^ 0) does not change this. To see this we may 
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perform the inverse in Eq. fl22|) to first order in M and then get an additional contribution 
of the terms 



-ika 



2i 



1 - 



-M + 



1 - 



+ 



1 - 



-iM- 



1 - 



iV 



-2iV 



-ika 



(36) 



to the S'-matrix. Now, both M and V have a factor of \fk~ c multiplying their cth columns 
(and rows) from their definitions and so a matrix element b cc i of the matrix in parentheses 
in Eq. (|3~6"D will have a \fk c and y/k^ dependence. An inelastic element of S(c ^ d) would 
now take the form 



S rr i 



-ik c a 



bed + 



T l/2 r l/2 
n Ac 1 Ac' 



-ikia 



E 



(r) 



E-iT x /2 i 



(37) 



As mentioned our interest is in the case when the entrance channel is at threshold so that 
this dependence is \fk~ e , making the inelastic probability |5 ce | 2 still continue to behave as 
k e ~ y/e. 



B. Overlapping Resonances 

Here we require the form of the S matrix near an energy E where many of the quasi- 
bound states may be simultaneously excited, i.e. the resonances overlap. Again, neglecting 
background for the moment, the S matrix is simply taken to be a sum over the various 
resonances. 

S = l-£ - Jf* - - (38) 



A 



E - E { { ] + iT x /2 



where A\ is a n x n rank 1 matrix with the cc'th component as T^T 1 /^ . There is no entirely 
direct justification of this form, but one can see that there is much which it gets correct. 

The A\ are symmetric, hence S is symmetric. Obviously it has the poles in the right 
places allowing the existence of decaying states with a purely outgoing wave at the resonant 
energies. A crucial additional assumption that also makes S approximately unitary is that 

1 /2 

the signs of the T Ac are random and uncorrelated both in the index A as well as c, regardless 
of how close the energy intervals involved may be. One simple consequence is that we 
approximately have that 

A X A X , = 5xx>T x A x (39) 

in the sense that the l.h.s. is negligible for A 7^ A' in comparison to the value for A = A'. 
With Eq. ( |39"D it is easy to verify the approximate unitarity of S. 

We investigate now the onset of the overlapping regime as E increases. D(E), the level 
spacing of the resonant E x , is a rapidly decreasing function of its argument. On the other 
hand, T x = T X1 + T X2 + • • • + T Xn , and since more channels are open at higher energy, T x is 
increasing with the energy of the resonance. The widths must therefore eventually overlap, 
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and T\ ^> D [E^j for the larger members of the sequence of E x r ' , s. In this regard there is 



a useful estimate due to Bohr and Wheeler [£L5[], that for n large 

T ) > ~ n . (40) 
D(EP) 

Appendix [A] derives this using a phase space argument. Here we point out that this is 
entirely consistent with the assumption of the random signs, indeed requiring it to be true. 
Take for example a typical inelastic amplitude 



-i 



r l/2pl/2 
1 Ac 1 Ac' 



a E^-E-iT x /2 



(c ^ d) (41) 



First let us note that the T\ being the sum of many random variables (the partial widths T\ c ) 
do not fluctuate much. Let T denote their typical value over the n overlapping resonances. 
Also since V = nD it follows that the typical size of a partial width T Xc is D. Therefore the 
typical size of the product r^ 2 r A {? is D but these random variables fluctuate randomly over 
the index A, and moreover the sign is random. Thus for energies in the overlapping domain 
S cc i is a sum of n complex numbers each of typical size D/T = 1/n, but random in sign. 
This makes for a sum of order 1 / ^fn. Clearly this is as required to make the n x n matrix 
S unitary. Note that the above argument fails (as is should) if c ^ d because then the signs 
of 

^Xc^Xc — Ta > are of course not random. 

Unlike the case of the isolated resonance, the S-matrix elements here are smoothly varying 
in E. Addition of a background term B cc i 

r l/2pl/2 

a E\' — E — iY x /2 

just shifts this smooth variation by a constant. If B cc i is also thought of as arising from a 
sum over the individual backgrounds then for the same reasons as discussed at the end of the 
preceding section |-B ce | 2 ~ k e ~ a/c for an entrance channel near threshold. For simplicity 
we will continue to take B cc r to be and look at the case with background in the appendix. 



VI. Q-MATRIX AND STICKING 

From the viewpoint of scattering theory, the sticking of the incident particle to the target 
is just a long-lived resonance. It is natural then to investigate the time-delay for the collision. 



Smith |T3] introduced the collision lifetime or Q-matrix 



Q = iUS— (43) 

which encapsulates such information. We review some of the relevant properties of Q. The 
rhs of Eq. fl43|) involves the 'open-open' upper left block of S so that Q is also an n x n 
energy-dependent matrix, having dimensions of time. For 1-dimensional elastic potential 
scattering S = and Q reduces to the familiar time delay ih 9 tj-^ . If v is a vector 
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whose entries are the coefficients of the incoming wave in each channel then ■u tr Q(£')'y is the 
average delay time experienced by such an incoming wave. Because physically the particle 
is incident on only one channel, v consists of all O's except for a 1 in the eth slot so that 
the relevant quantity is just the matrix element Q ee (E). Smith shows that this delay time 
is the surplus probability of being in a neighborhood of the target (measured relative to the 
probability if no target were present) divided by the flux arriving in channel e. This matches 
our intuition that when the delay time is long, there is a higher probability that the particle 
will be found near the target. 

Furthermore, as a Hermitian matrix, Q(E), can be resolved into its eigenstates v^ 1 ^ ■ ■ ■ 
with eigenvalues gi • • • q n . The components of are the incoming coefficients of a quasi- 
bound state with lifetime q\ and so on. Then 



n 
3=1 



(44) 



As can be seen from this expression, the average time delay results, in general, from the 
excitation of multiple quasi-stuck states each with its lifetime qj and probability of formation 
\v^^ -v\ 2 . However, we will find that using our resonant approximation to the S matrix near 
a resonant energy E^ the time delay will consist of only one term from the sum on the rhs 



of Eq. fl44j) , all the other eigenvalues being identically 0. 
Using equation Eq. ([43]). 



Q(J3)=ift E 



-iAy 



[E-EP-iTy/2 
which using Eq. (|39D simplifies to 

= E 



E 

AA' 



A X A X , 



h 



E 



(r) 



iTx/2 



E - 4? 



ir v /2 



(45) 



-.4 



A 



a (E - 4 r) ) 2 + 0V2) 2 

a remarkably simple answer. We need Q ee (E), where e is the entrance channel. 

fiTAe 

r(E- E^f + (r A /2)2 

r \(E-EPy + (T x /2y x rj 



(46) 



Qee(E) 



(47) 
(48) 



where the second equation has the interpretation (for each term) as the life-time of the mode, 
multiplied by the probability of its formation. Note how for each resonance E^ there is 
only one term corresponding to the decomposition of Eq. (0). The actual measured lifetime 
is the average of Q ee (E) averaged over the energy spectrum \g(E)\ 2 of the collision process. 



A. Energy averaging over spectrum 

With the target in state Q e (u) where c = e is the entrance channel, the energy of the 
target is fixed, and the time-dependent solution will look like 



12 



pi 00/ -ik c (E)x ik c (E)x 1 1 

i>(x, u,t)= dE\ g(E) ]T £ ce - -^=S(£) ce fi c («) . (49) 

7 \ c=i\JkjE) Jk^E) I I 



Recall, E is the total energy of the system. We are interested in the threshold situation 
where the incident kinetic energy of the incoming particle e — > 0. This can be arranged 
if g(E) is peaked at Eq with a spread AE such that i) Eq is barely above £ , * ar § et and ii) 
AE = 5e is some small fraction of e, the mean energy of the incoming particle. The second 
condition ensures that we may speak unambiguously of the incoming particle's mean energy. 
So, 



(Q ee (E)) = J dE\g(E)\ 2 Q ee (E) (50) 
dEQ ee (E) (51) 



1 

~ AE 

() denotes the average over the AE interval. Now Q ee (E) is just a sum of Lorentzians 
centred at the E^ s with width T\ and Eq. fl5"I|) is just a measure of their mean value over 
the AE interval. 

So long as the AE interval around which we are averaging, is broad enough to straddle 
many of these Lorentzians, the mean height is just 

x p(E)AE x (52) 



AE 



where the second factor is the number of Lorentzians in the AE interval and the third factor 
is the area under the 'Ath' Lorentzian. This is true regardless of whether or not they are 
overlapping. It will be convenient to write T\ as 

r A = nx2l A var( 7A ) (53) 

where var(7^) is the variance of the set of 7a c 's over the n open channels and k\ is a mean 
or effective wavenumber k c over the open channels, which for a particular realization A we 
take to be defined by Eq. (^3|) itself. Let ( ) denote the average over the occurrences of the 
quantity in the AE interval. T = k = (k\). Then Eq. (|52"D simplifies 



(Qee(E)) ~ 4 i^P^ ^) 

D nk{vax(^\)) 



ft < 55 > 



which tends to as k e ~ y/e. The form of Eq. ( |55| ) and all the steps leading up to it remain 
valid whether the Lorentzians are overlapping or not, as long as the AE = Ae interval which 
we are averaging over includes many of them. 

B. On an isolated resonance 

If the target is cold enough that the resonances are isolated, then as the incident particle's 
energy e — > 0, adhering to the condition Ae < e will eventually result in Ae becoming 
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narrower than the resonance widths. It becomes possible then for Ae to be centered right 
around a single isolated resonance at E { ^\ In this case {Q ee (E)) is found simply by putting 
E = E^\ because the spectrum \g(E)\ 2 is well approximated by S(E — E^). So 

«,_(£)) = ^ = = _-i . (56) 

Even in this case there is the y/e behavior as e — > and there is no sticking. 

In the extreme case that there are no other open channels at all (n = 1), (Q ee (E)) ~ 
r- because T\ = T\ e . In fact, e = 1, and (Q ee (E)) diverges, implying in this case 



that it is possible to have the particle stick. This is an exception to all the cases above but 
is experimentally not so relevant because we may always expect to find some exothermic 
channels open for a target with many degrees of freedom. 

VII. INELASTIC CROSS SECTIONS AND STICKING 

Another physically motivated measure of the sticking probability may be obtained by 
studying the total inelastic cross-section of the collision. The idea is that any long lived 
"sticking" is overwhelmingly likely to result in an inelastic colision process; i.e. that the 
scattering particle will leave in a different channel than it entered with. Using the original 
Wigner approach it is possible to show that for our case where we have only one scattering 
degree of freedom, the inelastic probability for an exothermic and endothermic collision 
vanishes like k e . The only possible exception to this is a measure zero chance of a resonance 
exactly at the threshold energy, El argct . In the event that there is a resonance E^ close to 
but above this threshold energy, it is only necessary that E is below E^ (by an energy of 
at least AE, the spread in energy) in order to observe the usual Wigner threshold behavior: 

^inelastic ~> like fc e OC yfl (57) 

for the inelastic probability. However our problem is unusual in the sense that because of 
the large number of degrees of freedom of the target, we will always find resonances between 
^target an( j j£ no ma ^ er ^qt^ sma ll E — _E* ar s ct = g is. Thus the Wigner regime is not 

accessible. Still the surprise is that a simple computation reveals the same behavior holds 
for large n: 

Pinclastic(£) = E P c-e(£) (58) 

= T,\Sce(E)\ 2 (59) 



= EEE Ac ( , e Ac r ; Ae — (60) 

c^e y 1 E E-^ iT\/2 E E^ + iT\/2 

— r 

=>- ^inelastic (E) = E n — ~ T Ae (61) 

v ' Y(i?-4 r) ) 2 + (r A /2) 2 1 ; 

where we used the random sign property of the T A { 2 's and the understanding that J2 Ta c — 
Tac = T\. Since the sum is over the n ^> 1 open channels, omission of a single term 

all c cj^e 
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can hardly matter. Apart from the factor h/T\, the rhs of the above equation is identical 
to the expression for Q ee (E) in Eq. fl48|). Averaging inelastic over many resonances E^ 
(overlapping or not) we may use the same algebraic simplifications as before to show 

(^inelastic) = J (62) 

As k e tends to 0, this gives the sfe Wigner behavior showing that there is no sticking. 

The above argument fails when there is only one open channel. There are no inelastic 
channels to speak of. In this case, if the energy E coincides with a resonant energy E-[' we 
will have the exceptional case of sticking, as discussed at the end of the previous section. 
But as pointed out there, this is primarily of theoretical interest only. 

VIII. CHANNEL DECOHERENCE 

The only case for which we stick is seen to be the case of when we are sitting right on 
top of a resonance with the incoming energy so well resolved that we are completely within 
the resonance width, AND there are no exothermic channels open. Having no such channels 
open amounts to an infintesimally low energy for a large target. Otherwise, the sticking 
probability tends to as y/e in every case. 

A. Time dependent picture 

From the time independent point of view, the physical reason for the absence of low 
energy sticking is contained in the factor ^ of Eq. ([48| ). This is the formation probability 
for the compound state. We will explain physically why it is small for n ^> 1. The resonance 
state is a many-body entangled state. If we imagine the decay of this compound state 
(already prepared by some other means say) each open channel carries away some fraction 
of the outgoing flux, with no preference for any one particular channel. Running this whole 
process in reverse it becomes evident that the optimum way to form the compound state 
is to have each channel carry an incoming flux with exactly the right amplitude and phase. 
This corresponds to however an entangled initial state. With all the incoming flux instead 
constrained to be in only one channel it becomes clear that we are not exciting the resonance 
in the optimal way and the buildup of amplitude inside is not so large; i.e., the compound 
state has a small probability of forming. 

The time dependent view is even more revealing. Imagine a wave packet incident on 
the system. For a single open channel Feshbach resonance, the build-up of amplitude in 
the interior region can be decomposed as follows. As the leading edge of the wavepacket 
approaches the region of attraction, most is turned away due to the quantum reflection 
phenomena. (It is a useful model to think of the quantum reflection as due to a barrier 
located some distance away from the interaction region.) The wavefunction in the interaction 
region constructively interferes with new amplitude entering the region. At the same time, 
the amplitude leaving the region is out of phase with the reflected wave, cancelling it and 
assisting more amplitude to enter. 
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Now suppose many channels are open. All the flux entering the interior must of course 
return, but it does so fragmented into all the other open channels. Only the fraction that 
makes it back into the entrance channel has the opportunity to interfere (constructively) with 
the rest of the entering wavepacket. The constructive interference is no longer efficient and is 
in fact almost negligible for n ^> 1, thereby ruining the delicate process that was responsible 
for the buildup of the wave function inside. The orthogonality of the other channels prevents 
interference in the scattering dimension. If we trace over the target coordinates, leaving only 
the scattering coordinate, most of the coherence and the constructive interference is lost, 
and no resonant buildup occurs. Therefore, one way to understand the non-sticking is to 
say that decoherence is to blame. 

B. Fabry-Perot and Measurement Analogy 

Suppose we have a resonant quantum mechanical Fabry-Perot cavity, where the particle 
has a high probability of being found in between the two reflecting barriers. Now, during the 
time it takes for the probability to build up in the interior, suppose we continually measure 
the position of the particle inside. In doing so we decohere the wave function and in fact 
never find it there at all. Alternatively, imagine simply tilting one barrier (mirror) to make it 
non-parallel to the first and redirecting the flux into an orthogonal direction, again spoiling 
the resonance. Measurement entangles other (orthogonal) degrees of freedom with the one of 
interest, resulting in flux being effectively re-directed into orthogonal states. Thus the states 
of the target (if potentially excitable) are in effect continually monitoring (measuring) to see 
if the incoming particle has made it in inside, ironically then preventing it from ever doing 
so. The buildup process of constructive interference in the interaction region, described in 
the preceding paragraph, is slower than linear in t. Therefore, the constant measurement 
of the particle's presence (and resultant prevention of sticking) is an example of the Zeno 
"paradox" in measurement theory. 

IX. CONCLUSION 

We have presented a general approach to the low energy sticking problem, in the form 
of i?-matrix theory. This theory is well suited for the task, since it highlights the essential 
features of multichannel scattering at low incident translational energy. We did not need 
to make a harmonic or other approximate assumptions about the solid target, which is 
characterized by its long range interaction with the incoming particle and its density of 
states. 'Warm" surfaces are included in the formalism, and do not change the non-sticking 
conclusion. 

Several supporting arguments for the non-sticking conclusion were given. Perhaps most 
valuable is the physical decoherence picture associated with the conclusion that there is no 
sticking in the zero translational energy limit. 

Reviewing the observations leading up to the non-sticking conclusion, we start with the 
near 100% sticking in the zero translational energy limit classically (sticking probability 
1). We then invoke the phenomenon of quantum reflection (Fig.|IJ), which keeps the incident 
particle far from the surface (sticking probability 0). Third, we note that quantum reflection 
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can be overcome by resonances (Fig. |2|), and since resonances are ubiquitous in a many body 
target, being the Feshbach states by which a partice could stick to the surface, perhaps 
sticking approaches 1 after all. Fourth, we suggest that decoherence (from the perspective 
of the incoming channel, with elestic scattering defmded as coherent) ruins the resonance 
effect, reinstating the quantum reflection as the determining effect. Finally, then, there 
is no sticking, and the short answer as to why is: quantum reflection and many channel 
decoherence. The ultrashort explanation is simply quantum reflection, but this is dangerous 
and non-rigorous, as we have tried to show. 

All this does not tell us much about how sticking turns on as incident translational 
energy is raised. This is the subject of the following paper, where a WKB analysis proves 
very useful. Quantum reflection is a physical phenomenon liked directly to the failure of the 
WKB approximation. 
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APPENDIX A: T ~ ND 

With the large number of degrees of freedom involved and assuming thorough phase 
space mixing associated with the resonance we may reasonably describe the compound state 
wavefunction by a classical ensemble of points (x,p x ,u,p u ) in the combined phase space of 
the joint system given by the normalized distribution 

i 

S(E - H(x,p x ,u,p u )). (Al) 



Pc(E) 



It is understood in the above that the system is restricted to be in the region x < a. This 
makes all accessible states of energy E with x < a equally likely. Then the rate of escape 
T/h through the hypersurface x = a of the members of this ensemble is 

-j— - [ dudp u [ dp x —5(E-H(x,p x ,u,p u )). (A2) 
i(E) J J m 



h p c (E) 



P3;g[0,Oo] 



p x /m is just the velocity in phase space of a point at x = a in the x direction. At x = a we 
have supposed no interaction. Hence the Hamiltonian separates in Eq. ([^). Therefore 



£ = 1 r ^ 7JpI\Jv ( p 2 

h p c (E) 



/^/n^jH^"(£ +Ftarset(M ' Pu) )) (A3) 

dudp u (A4) 



o 

1 



Pc 

JJtarget (y^p u )<E 

—n c ~ T^—tlg = -LnD. (A5) 
p c 2nhp Q * 2-kU 
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Therefore £ ~ n. pq (pc) is the quantum (classical) density of states (phase space volume) of 
the joint system at energy E. Qq (flc) is the quantum (classical) total number of states (total 
phase space volume) of only the target below energy E. We have used the correspondence 
between the Classical and Quantum density of states. 1/pq is identified with D, and the 
number of states of the target having energy less that E is just n, the number of open 
channels. 



APPENDIX B: INELASTIC PROBABILITY WITH BACKGROUND 

We show here that the inelastic probabilities remain essentially unaffected in magnitude 
with the presence of a background term in the S-matrix. In the isolated case the addition 
of b cc r to an inelastic element S cc ' simply changes the Lorentzian profile of jS^c'l 2 - In the 
more important overlapping case, the energy variation of S cc i is smooth in any case without 
background and 



I Srr' I 



r l/2 r l/2 
1 Ac 1 Ac' 



E 



(r) 



\B CC '\ 2 + X! 



E - iT x /2 
TaJac' 



E) +IY/4 



(Bl) 
(B2) 



where we have used the random sign property of the products T Xc i \, 



1/2 V 2 

r Ac to neglect the 2nd 
cross-term in comparison to the last one where again the same property is used to simplify 
the double sum to a single one. Summing over all the inelastic channels then leads to the 
same result of Eq. ( [IT]) with an added term of J2 \B CC >\ 2 which itself is proportional to k e 

Cy^e 



as discussed at the end of Section VB 
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